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Effect of Darling-Dennison and Fermi Resonance 
on Thermodynamic Functions 

Harold W. Woalley 

The effect on thermodynamic functions due to resonance of the Darling-Dennison type 
for water and of the Fermi type for carbon dioxide is calculated working directly with 
the elements of the secular determinants that define the energy levels. 

1. Introduction 

In calculating thermodynamic functions for various molecular systems, the complication 
is sometimes encountered that the energy levels are expressed indirectly by means of a deter- 
minantal equation. In many cases, the determinant will have nonzero elements only along 
and quite near to the principal diagonal. Frequently, no explicit general solution for the posi- 
tions of the levels as affected by the off-diagonal elements may be known, so that an explicit 
determination of the related thermodynamic properties as influenced by the effect in question 
would seem difficult at best. In some of these cases, useful results can be obtained in a general 
form from the secular equation without making use of explicit individual solutions. 

2. Secular Determinant Method 

If the secular determinant, \E ij —b ij W\, is expanded in polynomial form, beginning with 
the highest power of the unknown, W, then the successive coefficients for the different powers 
of Wgive the sum of roots and sums of products of the roots, such as ^W u 22WiW jf ZZ2WiWjW k , 

i i<j i<j<k 

etc., when allowance is made for alteration of signs in the successive terms. 

In the calculation of thermodynamic functions, the intermediate step of obtaining parti- 
tion functions involves the evaluation and addition of a large number of Boltzmann factors, 
using the energies of the various levels. This factor, e~ w i fkT , is expressible as a series so that 
the summation over a set of energy levels 

S e-^'^S ( ~ k Z )m S W? (1) 

may be found if the individual sums of powers of energy 2W*, 2"PP?, 2 Wf, etc., are known. It 
is possible to find these quantities in terms of sums of products of the roots and hence in 
terms of the coefficients of the polynomial expansion of the secular determinant. In this process 
several algebraic steps are involved including the application of the multinomial theorem. 
The result can thus be put in the form 

S w?=GC Wif-2 ss WtW, (2) 

i i i<j 

S w?=(S w t y-z CS iiiXSS w,w k )+z sss w t w,w t (3) 

i i i j<k i<j<k 

2 w?=(E M',.)<-4(S TWXJS W,w*)+4: (2 ^)ESE WjW*w,) 

i i i j<k i j<k<l 

+2(Si; w t w,f- 4SSSS^^^ w, (4) 

i<i i<j<k<l 

i i i j<k i j<k<l 

+5 (S ^)(SS WjW k y-5 (S ^CSSS w,w t w,wj 

i j<k i j<k<l<m 

-5 es w t Wj) ESS WiW.^+sssSSS w t w i w k w l w m (5) 

i<j k<l<m i<j<k<l<m 
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i i j < fc i j < k < I 

i j < k i j < /c < I <w 

- i2(s wo(ss wwoess h^.^. w 7 ») 

i j < fc Z < m < w 

« i < ft < J < m<n 

i < j i < j k < I < to < rc 

+3(ESS^^^) 2 -6I]SSSS2^^^F i if w ^). (6) 

i <j<k i <j <k < l<m<n 

These coefficients for the various terms of the polynomial expansion, y^.y^. WtWu 

i <j 

SSS^^Aj etc > can be obtained in the case of a specific secular determinant and 

i <j <k 

introduced into eq (2) to (6) . If the secular equation 

{Eti-dtsWl^O (7) 

involves a determinant of the nth order with nonzero elements along the principal diagonal and 
along diagonals immediately adjacent to it, then it is possible to represent the coefficients 
explicitly and introduce them into eq (2) to (6) so as to represent the quantities 2WP directly 
in terms of the elements of the determinant. Thus for the symmetric determinant with E u = 
k i} Ei t i + i=E i+ i,i=<ii, E if i +2 =E i+ 2,i = b i , E iti ^ z =E i+z ,i : =Ci J and with all other elements equal to 
zero, one may with considerable algebra transform from the previous equations to the following: 

Z)^ = Z)*< (8) 

i i 

SW^=Sfc?+2&tf+ 2 &?+2iSc? (9) 

i i=l i = l i'=l i=l 

n 7i—l n—2 71—3 

i i=l imLk* i=1 * = l i=1 

+ 6 i E,a i a i+1 b i + QT l a i b i+1 c i + 6J2a i+2 b i c i , (10) 
1=1 <=i 1=1 

including terms involving a's, b's and c's. If the c's are zero, then 

i i=l i=l i=l 

i=l i=l i=l i=l i=] 

i = l i=l i=l i=l t = l 

and if the 6's also are zero 

SWJ=S*f+BSaf(*?+*f**+i+***f + i+tf + i) 

i i 

+ 5£V 4 (^ + ^ + i) + 5Safa? + i(^ + 2/: i+1 + ^ +2 ) (12) 
i i 
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and 

ill 1 

+6S« +1 (i?+2t l i i+1 +i i i <+2 +3t? +1 +2i i+1 i <+2 +i? +2 ) 

i 

+ 2Sa?+6Sa?af +1 + 6Sa?o?+i + 6Sa?a? +1 a? +2 . (13) 

ii i i 

Equations (8) to (13) may also be obtained by a matrix method with less effort. This is 
based on a theorem of Borchardt [l] 1 and of Sylvester [2] which states that if A is a matrix then 
the characteristic roots of A n are the 7ith powers of the characteristic roots of A. Then if H is 
the energy matrix, the sum of nth powers of roots W t of the secular equation is given by the 
trace of the nth power of the matrix H: 

Y}W?=tmceH n . 



(14) 
(15) 



so that eq (1) may also be written as 

£>-wv* t = trace e -ff/tr 

i 

and, in effect, include all the results of eq (8) to (13). 

3. Application of the Method 

Among cases in which molecular energy is given by a secular equation are resonances of the 
Darling-Dennison type and the Fermi type. In the case found by Darling and Dennison [3] 
for the H 2 molecule, occurring also in some other cases, the energy W is given by the secular 
equation 

Wl-W W 12 



W„ 






wi-w W 23 

w 32 Wl-W w 3l 







w„ 



wi- w 



= 0, 



(16) 



where the energies for the individual interacting levels would have been W° if the perturbation 
had been negligible. The magnitude of the perturbation is indicated by the off-diagonal terms 



^,m=WV 1 , i =C^^3+2=i7V^(^-l)(^3+l)(^3+2). 



(17) 



Any state indicated by v u v 2) v* with #i>l is perturbed appreciably by a state ^ — 2, v 2 , #3+2, 
because of the approximate equality of v x and v%. For higher quantum numbers, large groups 
of levels will be perturbing each other, and the secular equation will involve a secular deter- 
minant of correspondingly high order as indicated by eq (16). In the discussion which follows, 
the number of levels perturbing each other in a single group will be indicated by n, which will 
also be the order of the determinant. The mean value of the energy for the group of levels 
expressed as wave numbers will be indicated as G. The increment in the vibrational partition 
function due to the presence of the perturbation for a group of interacting levels is given by 

AQ = e -> cOI *r [_£ A?TFi+ l Q0 A?w? _ 1 Q0 A ^ w? . . J (lg) 



'.Figures in brackets indicate the literature references at the end of this paper. 
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According to the formulas derived and the definition of G, the various increments in the sums 
of powers of the energy W measured from the mean value will be given by 



A j2Wi = \ 



A2FN2S^, i+1 \ (19) 

A£W?=32^, i+1 (^+C +1 ), etc. I 

i i / 

In the remainder of the treatment for this case of resonance, only the effect of A^ W° on 

i 

the partition function and derived quantities will be considered. There is a basis for regarding 
higher order effects, as of A^JFf, etc., as negligible. The levels which are involved in a given 

i 

perturbing set have the same value of Vi~\-Vz and are subject to the additional restriction that 
either Vi or #3 is only odd or only even within a single group. There are thus four cases to con- 
sider: 

I. Vi ranges from 2n — 2 to while Vz ranges from to 2n— 2. 

II. Vi ranges from 2n— 1 to 1 while Vz ranges from to 2n— 2. 

III. Vi ranges from 2n— 2 to while Vz ranges from 1 to 2n — 1. 

IV. Vi ranges from 2n—l to 1 while Vz ranges from 1 to 2n— 1. 
These four cases give, respectively, for W? ii+1 

I. Wl i + 1 =^(2n-2i)(2n-2i-l)(2i-l)(2i) (20) 

II. Wli +1 =j (2n+ 1 - 2%) (2n-2i) (2i- 1) (2i) (2 1) 

III. Wl i+1 =^(2n-2i)(2n-2i-l)(2i)(2i + l) (22) 

IV. W 2 hi+1 =^ (2n+l — 2i) (2n-2i) (2i) (2i+ 1), (23) 

n-l 

where in each case i ranges from 1 to n— 1. The evaluation T^, W* ti+1 is quite straightforward 
and gives for these cases, respectively, 

y 2 2n(n— l)(n+l) (A 2 . , AS 

~ zjr (4n — IOti+9) (24) 

„2 



^ 2n(n-l)(n+l) (2 ^ 1)(2n+1) (2S) 

7 2 2ti(ti-1)(?i+1) 



4 15 

and 

7 2 2n(yi-l)(7i+l) 
4*~ 15 



(2ti-1)(2ti+1) (26) 

(4n 2 +10n + 9). (27) 



Considering the molecular vibrational energy as if it were due to harmonic oscillators 
having frequencies «i and co 3 wave numbers, for the two vibrations it is found that } the average 
energy for the group of n levels is given in the various cases by 

I. 0=(n-l) («H-««) (28) 

II. G=n^ + (n- 1)«3 (29) 
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III. G=(n-l)<* x +nm (30) 

IV. G=n(u x + u 2 ) (31) 

The expressions are now at hand to permit the estimate of the relative increase in the 
partition function due to the perturbation. Continuing to treat the various cases separately, 
there are obtained for case I, 

V ), Ac he 

03 = 01=0 

,^(j^Y d -^ ai X 1 -«"? ") [,-ft **«> + 6^ ( ^ ) +9e-^ ( »^1 (32) 
V 7 [ 1 _ e -^ ( "» + " 3, J 

for case II, 

f£«a^ -7? <-i+»»)» / ftc V y 3 2n{n 2 - 1) (4rc 2 - 1) 
m\ J hi" \kT) 4 15 

03 = 01 = 

( hc Al-e-£ ai )(l-e-^)e^ l -2^KW -3^«o 1+ . 3 ) ^ch-,)] 

(33) 
for case III, 

TJft v -£f <»i+»»)« / frc V 7l 2w(w 2 - l)(4n 2 - 1) 



( 



m?\ = e kT S e u?v 4 " is 

O )o he he 

01 = 03 = 

v y [ 1 _ e -ii ( " i+ " 3) ] 

(34) 
and for case IV, 

s ^>.,., ( | ) -^ wi „ t „ 



veA 



'1=0 »! = 



= 7 2 (ir^) V1 ^ -, 6 ; l_9e * r +6e * r ' 3 +e kT ' 3 'J- (35) 

V y [l_ e -ll ( " 1+ " 3) ] 
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Then, in lowest order of approximation, the total effect of the perturbation is indicated by the 
sum of the separate (AQ/Q)t for the four cases, to give as the approximate value of A In Q 



A\nQ = ^(AQ/Q) l = y i Qf\ 2 ^ 



hc^ 

' kT 



-)(.- 



he 

"Fr* 



A 



~um ( w l+ w 3) 



- I g"^ (a,1+ " 3) 



„ he , . „ he , , . he . . . 



+ Ve 






■ + « 



H( 



? -^K^3) +10e -2^(^+. 



> + 3e -^^ 



)l 



(36) 



If the approximate equality of a>i and co 3 is now taken into account, then letting u=hco)/kT 
represent either hcui/kT or hcuz/kT or their average value, it is found that 



AlnQ = 



y 2 u 2 (l-e~ u ) 



(1- 



2 ^[e- 2u + Qe- 3u + 15e-' u + 20e- 5u + 15e- Qu + 6e- 7u + e- 8u ] 



7 2 u 2 e~ 2u 



'a> 2 (l-e- u Y 



(37) 



This is the approximate contribution to the Gibbs free-energy function, —F°/RT, due to 

the effect of the perturbation. The corresponding contributions to H°/RT and C°/R are 

d d 2 

—u -j- A In Q and u 2 -j— 2 A In Q, which give 



and 



^/M0^tf & ^.2^g±^-l} 



A(C°JR)=^u 2 



(l-e~ u y 



2 U-4u 



d-6-") 



-2w 



(l + 3e- u + e- 2u ) 
(l-e- u ) 2 



(38) 



(39) 



A value of 7 of 77.52 cm -1 has recently been given by Benedict [4] with higher order constants 
and has been used in the present calculation. For the present case, o> may be taken as about 
3,650 cm -1 . As should be expected, the total effect is found to be small. This is due largely 
to the fact that the lowest levels showing this type of resonance are at about 7,300 cm -1 so 
that a relatively high temperature is required before the levels in question are appreciably 
excited. The magnitude of the effect is indicated in table 1 for several temperatures. 

One may inquire as to the approximate behavior of the effect in the limiting condition of 
high temperatures. If one neglects higher order effects in the manner of the foregoing deriva- 
tion, which would neglect an appreciable relative effect due to anharmonicity, then when the 
functions are expanded in increasing powers of u, the leading term gives the approximate values 



- AF°/RT « -4-w AH°IRT « 2 
artr 



and AC°JR=6 



OfW 



(40) 



or 



Table 1. Effect of Darling-Dennison resonance on thermodynamic functions for H2O 



Temperature 



1,000. 
2,000 _ 
3,000. 
4,000. 
5,000- 



A(—F°/RT) 



2.6X10-7 
1.88X10-5 
0.784 X10-* 
1. 75X10-* 
3.03 X10- 4 



A(H°/RT) 



24X10-7 

7.81X10-5 

2.37X10-* 

4.52X10-* 

7.23X10"* 



AC$/R 



206X10-7 
33.1X10-5 
8.01 X 10-* 
14.2X10-* 
22. IX 10-* 



-AF'/RT-t^AHVRT^^ and A<7;/B~6 £Q£f- (41) 
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The quadratic power of temperature is a result of the fourth power dependence of AS W 2 upon 
the quantum numbers, which supplies a fourth power of T as a factor in the summation. This 
is cut down to the second power by the inverse second power in the coefficient of this term. 

Another case of resonating levels for which details have been worked out is that of Fermi 
resonance in the case of C0 2 . While the effect of the perturbation upon the thermodynamic 
functions can be expected to be fairly small in this case also, since the perturbation, as in the 
previous case, involves the displacement of some levels upward and some downward in such a 
way as to nearly compensate, it is nevertheless a larger effect than in the case of H 2 0, mainly 
because the energy of the lowest resonating levels is near 1330 cm -1 for C0 2 as compared with 
the 7,300 cm" 1 for H 2 0, so that comparable excitation takes place at a lower temperature. 

It was shown by Dennison [5] and Adel and Dennison [6] that a satisfactory representation 
of the C0 2 spectrum could be obtained if account were taken of perturbations within groups of 
levels with equal values of I, the azimuthal quantum number, and with equal values of 2vi J {-v 2 . 
Then the effect of the perturbation within a group of n levels of this kind is indicated by the 
shifts of roots of such a secular equation as (16) with corresponding effects on the thermodynamic 
functions indicated in terms of the partition function according to eq (18) and (19). In each 
group of n levels, the value of v { ranges from zero to n— 1 , with v 2 going from its maximum value 
to its minimum value which is equal to |/|. The number of levels interacting within one group 
is given by n=^(V— \l\) — 1. The elements of the determinant which are adjacent to the main 
diagonal were presented by Dennison in a form which can be so converted as to indicate that 



Wl i+l =W 2 i(n-l)(V+2-n-i), 



(42) 



where i ranges from 1 to n— 1. For the evaluation of the principal contribution to the parti- 
tion function according to eq (18), there is needed a knowledge of the distribution of the 
unperturbed vibrational levels. It is assumed that these are given in the somewhat uncon- 
ventional form 

ft-E^^+2^^:fe-i)+S^^ (43) 

Then the average value of G for a given group of levels is given by 

G=o: 2 V+co s v 3 +c, n (n-l)+X 22 (V 2 -V)+X 3 M-v s )+X nn (n-l) 2 

+X 23 Vv 3 +X Vn V(n-l)+X n3 (n-l)v 3 , (44) 



where 



o) n =—o)i — co 2 — ~X n + — X 22 — -Xi 2 



X nn — —X\\-\-—X 22 — ~^.X\ 2 



X vn — — Xn — 2X 22 



X n3 —^X n - 



-x 2: 



(45) 



By simple summation, 

and for the next term of (18) 

r . 2 ft(n 2 -l)(4-n : 



T,wi i+l m+wuj=wi- 

! = 1 



30 



[a„-8,V-X na v& 



(46) 



(47) 
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12 1 

where 8 V =— X n — - ( zX 2 2— ^ ^12- After multiplying by the exponential Boltzmann factor for the 

average vibrational energy, the summation over V and n follow, but with the requirement that 
half values of the summand are to be used when Fhas its lowest value, 2{n — 1), in order to take 
account of the absence of I degeneracy for £=0. The choice of doing this rather than multiply- 
ing by 2 for 1 5*0 comes about as the result of the fact that half of the rotational levels are 
absent due to the symmetry of the molecule and the zero nuclear spin for oxygen. Then this 
is compared with a total approximate Q similarly omitting any factor for / doubling. If the 
effects for anharmonicity are disregarded, the leading term is obtained as 



AQ l/hc\ a WU 



2hco} 2 / 

k~ 



C«2 / hC(>)i\ 



Q 2\kTJ / _2ftco L2 y/ ftcu, 2 \: 

Vl-e kT )\l-e kT ) 



(48) 



The effect of Fermi resonance, using u=hoa 2 /kT and disregarding the distinction between w 1 
and 2<o 2 , is thus given approximately [7] by 



A(-F°IRT) = 



\Wl 



2 co* (l-e- 2 ")(f-e-") 2 



and 



C%iB=A(-F°IRT).2\2[u^ 



2e 



(l-e- u ) 2 ' (1-e 



- 2 ") 2 J 1 \' 



(49) 
(50) 

(51) 



The effect can be obtained in more detail by including an additional contribution from eq 
(46). Using Ui=hco)i/kT and using u x = 2o) 2 after the first term, eq (48) is replaced by 

A{nq ~2\kTJ (l-e- 2u 2) 2 (l-e- u 2) 2 

"QfJ ^{ 2(1 -e-***)* [e- 2 ^+2e-^ + Qe-^ + Qe~^+Se-^] 
X t 



■2u^ [e~ 2U2 + 2e- 3u *+ 1 7e-* u 2+20e- 5u * + 33e- Qu 2+ 1 8e~ 7u 2+ 9e- Su *] 



[ e -2u 2 j_ 4^-3^2+19^-^2 + 3 0e~ 5u 2 + 396- 6w 2+246- 7w 2 + 96~ 8w 2] 
[ 6 -2w 2 + 4 6 -3w 2 + Se~ iu 2 + 8e- 5w 2 + 36~ 6w 2] 



(l-e~ 2w 2) 

X 23 e~ u * 



(l-e- u 3)(l-e- 2 «2) 
X n3 e- U * 



2(1— 6~ M 3)(i— e~ 2u 2) 
-X22 



_J 6 -2M 2 +2^- 3W 2+6e- 4W 2+66- 5M 2+36- 6W 2] 



(1— e - 2u *y 

B+g 22 



6(l-e- 2w 2) 



[e- 2u 2+l0e~ 3u 2 + 4:3e-* u 2 + 88e- 5u 2+103e- Qu 2+70e- 7u 2 + 21e- Su 2] 

[e-3t*2+8e- 4tt 2 + 4l6- 5tt 2+i04€- 6tt 2+1196- 7M 2+926- 8tt 2+i9g- 9tt 2] 

3X 22 e~ 4u 2 X 12 e~ 5u 2 



{l-e- 2u 2)\l-e- u 2f (\- e - 2u 2) (1 - e - u if (1 - e - 2w 2)2(i — e -u 2 y 
X l3 e- U 3e~* u 2 X 23 e- U 3e~ 3u 2 



2(1 -6~ M 3) (1 _g-2« 2 )2(l — 6 -«2)2 (1 - e -U 3 ) (1 _g-2« 2 ) ( X 



— <?-»2)3 ) 



(52) 
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A calculation was made for C0 2 using as values for the constants (in cm" 1 ) 
«!= 1342.86 co 2 = 666.27 « 8 =2349.15 
A^ n =-2.20 X 22 =-0.75 Xm=- 12.46 



X 2 =3.76 



X 13 = - 18.54 X 23 = - 12.48 



with B+g 22 = 1.42 and W =51.01 cm" 1 , giving a>„=3.39, X n 



-4.24, X^=3.38, and Z n3 =3.21. 



With these constants, the magnitude of the effect as given by eq (48) is indicated in table 2 as 
the first column under each heading while the combined effect consistent with eq (52) is shown 
as the second column under each of the headings, A(—F°/RT), A(H°/RT) and AC° P /R. The 
entire effect for C0 2 is rather small and in many computations could be neglected. For elevated 
temperature, the leading term in this case is roughly linear in the absolute temperature rather 
than quadratic as was found to be the case for H 2 0. This is due to the essentially cubic 
dependence of Wf ti ^ 1 upon quantum numbers in the present case, as compared with the 
quartic dependence in the case of H 2 0. It is also seen that this simple dependence upon tem- 
perature is modified greatly by the effect of anharmonicity. The large effect of anharmonicity 
at high temperature can be regarded as due to its lowering of entire groups of resonating levels 
with a resultant increased population in such groups, which provides an increase in thermody- 
namic effect. 

This treatment does not make any adjustment for the effect of resonance on the composi- 
tions of states, since the rotational perturbation is quite small, and, in addition, the resonance 
shifts have compensating effects in first approximation as has been noted in deriving the larger 
contributions due to perturbation which have been taken into account. They would thus be 
quite negligible. 



Table 2. Effect of Fermi resonance on thermodynamic functions for C0 2 



Temperature 

°K 

200 

250 

300 

400 

500 

1,000 

2,000 

3,000 

4,000 

5,000 



A(-F°/RT) 



A(H°/RT) 



MCp/R) 



5X10-6 
22X10-6 
5X10-5 
18X10-5 
33X10-5 



36X10-6 
12X10-5 
26X10-5 
56X10-5 
85X10-5 



27X10-5 
70X10-5 
11X10-* 
17X10-< 
22X10-< 



eq(49) 



0.0012 
.0029 
. 0045 
.0061 
.0077 



eq(52) 



eq(50) 



0. 0012 
.0030 
.0048 
.0067 
.0087 



0. 0018 
. 0033 
. 0048 
.0063 
.0078 



eq(52) 



eq(51) 



0. 0018 
. 0036 
. 0056 
.0078 
.0103 



0. 0033 
. 0062 
. 0093 
.0124 
.0155 



eq(52) 



0. 0035 
. 0073 
.0119 
.0173 
.0233 



Equations (49), (50), and (51) also represent approximately the effect of Fermi resonance 
in H 2 involving levels for which 2z?i+?; 2 is a constant, even though v 2 is not degenerate as for 
C0 2 . Benedict [4] gives the matrix element (v u » 2 +2|»i + l, vd = %{200(vi + l)(v 2 + 2)(v 2 + l)}*, 
indicating that Wl— 200 cm -2 . A few values for the Fermi resonance effect for H 2 are given 
in table 3. The effect is much smaller than the Darliiig-Dennison resonance effect at the higher 



Table 3. Effect of Fermi resonance on thermodynamic functions for H 2 



Temperature 

°K 

500 

1,000 

2,000 

3,000 

4,000 

5,000 



A(-F°/RT) 



A(H°/RT) 



8X10-8 


10-8 


3X10-6 


8X10-6 


10-5 


2X10-5 


2X10-5 


3X10-5 


3X10-5 


4X10-5 


4X10-5 


4X10-5 



ACIIR 



4X10-8 
2X10-5 
4X10-5 
5X10-5 
7X10-5 
9X10-5 
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temperatures although it is the larger at 1,000° K and below. Its earlier rise is partly due to 
the resonance occurring first near 3,400 cm" 1 as compared with 7,300 cm -1 for the Darling- 
Dennison resonance. The magnitude of the resonance effects for H 2 and C0 2 are shown in 
figures 1 and 2, respectively. 



.002 
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1 1 




- 


V 


.001 
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-^^^£hX 


' 1 T 





.0 2 




1 1 1 


1 
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- 




1^^ 







""l 1 1 


1 



12 3 4 5 

T/lOOO 

Figure 1. Effect of resonance on thermodynamic 
functions for H 2 0. 

ci, hi, fi: Contribution of Darling-Dennison resonance to C v jR, 
(H°-E§)/RT, and -{F°-E%)IRT. c, h, f: Contribution of 
Fermi resonance to C$/R, (H°-E )/RT, and -(F°-E%)/RT. 



1 2 3 4 5 

T/lOOO 

Figure 2. Effect of Fermi resonance on thermody- 
namic functions for C0 2 . 

c,h, f: Contributions to C./R, (H°-E° )/RT, and -(F°-E° )jRT. 
Upper branches with anharmonicity, lower branches neglecting 
anharmonicity. 



4. Conclusion 

A useful method of calculation of thermodynamic effects for levels denned by a secular 
determinant equation has been demonstrated. 

Many other instances with similar perturbations are known and their effects on thermo- 
dynamic functions should be calculable by similar methods if the theoretical description of the 
perturbation has been worked out and the spectra fitted in the detail that has been achieved 
for H 2 and C0 2 . In particular, C 2 H 2 gives an example of Darling-Dennison resonance for 
which eq (37), (38), and (39) would apply. Other molecules with Fermi resonance similar to 
that for C0 2 to which eq (49), (50) and (51) apply include CS 2 , C 2 H 6 , C 2 D 6 and N 2 0. 
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